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Abstract 

Electronic coherence dynamics in irans-polyacetylene oligomers are considered by explicitly com- 
puting the time dependent molecular polarization from the coupled dynamics of electronic and vi- 
brational degrees of freedom in a mean-field mixed quantum-classical approximation. The oligomers 
are described by the Su-Schrieffer-Heeger Hamiltonian and the effect of decoherence is incorporated 
by propagating an ensemble of quantum-classical trajectories with initial conditions obtained by 
sampling the Wigner distribution of the nuclear degrees of freedom. The electronic coherence of 
superpositions between the ground and excited and between pairs of excited states is examined 
for chains of different length, and the dynamics is discussed in terms of the nuclear overlap func- 
tion that appears in the off-diagonal elements of the electronic reduced density matrix. For long 
oligomers the loss of coherence occurs in tens of femtoseconds. This timescale is determined by 
the decay of population into other electronic states through vibronic interactions, and is relatively 
insensitive to the type and class of superposition considered. By contrast, for smaller oligomers the 
decoherence timescale depends strongly on the initially selected superposition, with superpositions 
that can decay as fast as 50 fs and as slow as 250 fs. The long-lived superpositions are such that 
little population is transferred to other electronic states and for which the vibronic dynamics is 
relatively harmonic. 
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I. INTRODUCTION 



Electronic decoherence (the decay of the off-diagonal elements of the electronic reduced 
density matrix) in molecules is a basic feature of the electron- vibrational evolution that ac- 
companies photoexcitation [1], passage through conical intersections [2], energy transfer {^J 
or any other dynamical process that creates electronic superposition states. In the decoher- 
ence language 4J, the electrons are the system of interest, the nuclei act as the bath, and it 
is the system-bath interaction what leads to the decoherence. Establishing mechanisms for 
electronic decoherence is central to our understanding of the dynamics underlying fundamen- 
tal processes such as photosynthesis, vision or electron transport . It is also vital in 



a, 



the development of approximation schemes to the full vibronic evolution of molecules 
and it is the starting point for the design of methods to preserve the coherence of electronic 
superpositions in molecules that can be subsequently exploited in intriguing and potentially 
useful ways via quantum control [8] or quantum information Q schemes. 

Timescales for electronic decoherence in polyatomic molecules are often exceedingly fast, 



on the order of tens of femtoseconds I1|.|10H12||. This timescale is normally determined by the 
vibrational degrees of freedom of the nuclear dynamics, with slower torsional, rotational or 
possible solvent dynamics (if present) playing a secondary role 10|. However, challenges in 
understanding electronic decoherence have arisen from recent spectroscopic observations that 
have demonstrated that in some photosynthetic systems electronic coherences can be long- 



lived 
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15] , with lifetimes exceeding 400-600 fs. These results have lead to discussions of the 



role of quantum coherences in biological processes and reconsideration of our understanding; 

24]. Many 



25| 



of decoherence dynamics in single molecules and molecular aggregates 0, [13L S 
of the associated computations utilize phenomenological models or master equations 
that approximate the dynamical effects of the bath on the system coordinates without 
explicitly following the bath dynamics. In these approaches, the effect of the bath on the 
dynamics is typically determined by adjustable parameters that can be chosen to reproduce 
experimental findings, when available. When possible, however, explicitly following the 



dynamics of the nuclei is much preferred 26j, |27|. This is because electronic decoherence in 



molecules can be understood as arising from nuclear dynamics on several electronic potential 
energy surfaces fll. [lol 12 1. For example, for an entangled vibronic state of two levels, of the 
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form 

\m) = Mxim+\<i>j)\xj(t)) a) 

where the \(f> n ) are orthonormal electronic states and \Xn{t)) denotes the nuclear state in the 
n th electronic surface, the electronic reduced density matrix p e is given by: 

Pe (t) = Tr N {Mt))m)\} 

(2) 

= l&)(0i|<xi(*)lxi(t)> + IhXhKxjWMt)) + [IfoKhKxMxiit)) +h.c] . 

Here the trace is over the nuclear states and h.c. denotes the hermitian conjugate. Hence, 
the decay of the off-diagonal matrix elements in p e (t), i.e., electronic decoherence, is governed 
by the degree of overlap of the nuclear wavepackets = (Xj(t)\Xi{t)) associated with 

the electronic states in the superposition. Thus, by understanding the events that lead to 
a decay of the overlaps (t) one obtains direct insights into the mechanism of electronic 
decoherence between states i and j. A schematic representation of such evolution and decay 
for a particular pair of states is shown in Fig. [TJ 



E 
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FIG. 1: Evolution and decay of the overlap of the nuclear wavefunctions in two electronic potential 
energy surfaces Sij(t) = {Xj(,t)\Xi(t)) upon instantaneous excitation from state \4>i) to state \4>j). 
Both anharmonicities and population transfer to other electronic states (dotted lines) can lead to 
a decay of Sij(t) and thus to decoherences between \<j>i) and \(f>j). In the scheme, E is the energy 
and q denotes the general nuclear conformational space. 

In this paper we present a study of the electronic coherence dynamics in trans- 
polyacetylene (PA) oligomers in which the dynamics of both electronic and vibrational 
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degrees of freedom are explicitly taken into account. We do so in an approximate scheme 
where the nuclei are considered classically and the electrons quantum mechanically. The 
oligomers are described using the well-known Su-Schrieffer-Heeger (SSH) Hamiltonian 28 ]. 
The SSH model treats the molecule as a tight-binding chain in which the electrons are cou- 
pled to distortions in the polymer backbone by electron-vibrational interactions. In spite 
of its simplicity, the SSH Hamiltonian is remarkably successful in capturing the electronic 
structure of PA, its photoinduced vibronic dynamics and the rich photophysics of polarons, 
breathers and kinks 29^ 31] . This model is often used to study the dynamical features 



caused by strong electron-ion couplings 



33]. 



The coupled dynamics of nuclear and electronic degrees of freedom of the molecule is 
followed in a mean-field (Ehrenfest) mixed quantum-classical approximation 34j-l36| and de- 
coherence effects are incorporated by propagating an ensemble of quantum-classical trajecto- 

tion AH-, 



ries with initial conditions selected from the nuclear Wigner distribution function 



39] 



of the chain. In this way the dynamics reflects the initial nuclear quantum distribution 
and is subject to the level broadening and internal relaxation mechanism induced by the 
vibronic couplings. Using this model we study the possible effect of system size, nuclear 
initial conditions and type of electronic superposition states on the dynamics of electronic 
coherence. 



II. MODEL AND METHODS 



A. The SSH Hamiltonian 



The SSH Hamiltonian |28[ models PA oligomers as one-dimensional tight-binding chains, 
each site representing a CH unit. The Hamiltonian for an A-membered oligomer is given 
by: 

-ffssH = H e i ec + H ph , (3) 

where 

JV-l 

H e l ec = 2j ^2 + _ U n)}(4i+l,s C n,s + 4i,s C n+l,s) and, 

n=l s=±l (4) 

?1=1 71=1 
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are, respectively, the electronic and nuclear parts of the Hamiltonian. Here, u n denotes 
the displacement of the nth CH site from the perfectly periodic position x = na with a as 
the lattice constant, M is the mass of the CH group, p n is the momentum conjugate to u n 
and K is an effective spring constant. The operator (or c n>s ) creates (or annihilates) a 
fermion on site n with spin s and satisfies the usual fermionic anticommutation relations. 
The electronic component of the Hamiltonian consists of a term describing the hopping of 7r 
electrons along the chain with hopping integral to and an electron-ion interaction term with 
coupling constant a. The quantity a couples the electronic states to the molecular geometry 
and constitutes a first-order correction to the lowest-order hopping integral to- Throughout 
this work, we use the standard set of SSH parameters for PA: t = 2.5 eV, a = 4.1 eV/A, 
K = 21 eV/A 2 , M = 1349.14 eV fs 2 /A 2 , and a = 1.22 A. 



B. Ehrenfest electron-vibrational dynamics 



The electron-vibrational dynamics of^t 
mixed quantum-classical approximation 



he chain is followed in the mean-field Ehrenfest 
jl. [35L 38]. In this approximation, the nuclei move 



classically on a mean-field potential energy surface with forces given by 

p n = -(m^ L Mt)}- (5) 

In turn, the antisymmetrized M electron wavefunction \<p(t)) satisfies the time-dependent 
Schrodinger equation 

ih^\<p(t)) = H SSK [u(t)Mt)), (6) 

where u = (m, U2, ■ • • , un)- Since Hssu is a single-particle operator, the electronic properties 
of the system are completely characterized by the single-particle electronic density matrix 



A 



^t) = ^(^\ c lsC m ,Mt)). (7) 



From Eq. ((6]) it follows that the dynamics of p n ^ m satisfies 

ih^p n! m{t) = y^(<l(*) I [4,s C m,s, H elcc \\Lp(t)) = y^X h rn,m'Pn,m'(t) - /i m >p m ', m (*) ) , (8) 
s m' 

where h n>m = (n, s\H e i ec \m, s) are the single-particle matrix elements of H e \ ec and \n,s) = 
c]j S |0) where |0) is the vacuum state. 
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Equation (jSj) is integrated by decomposing p n , m (t) into orbitals. For this, let \e, s) be the 
eigenorbitals of spin s and energy e of the system at preparation time (H e \ ec (t = 0)|e, s) = 
e\e, s)). Using this basis, the initial electronic reduced density matrix can be expressed as 

N 

Pn .in 

e,e'=l s 

where («/p(0) |c| s c e > iS \ip(0)) characterizes the initial electronic distribution among the single 
particle states, and \e,s) = c| s |0). In writing Eq. Q we have employed the basis transfor- 
mation function s = J2 e= x(c, s\n, s)c\ s . We adopt the ansatz that upon time evolution 
p n ,m(t) maintains the form in Eq. ([9]). That is, 

N 

Pn, m (t) = S ' n ' S ^ m ' S ' e 'W' S )(^(°)l<^',^( ))- ( 10 ) 

e,e'=l s 

The utility of this ansatz is that if the time-dependent orbitals |e(t) , s) satisfy the single- 
particle Schrodinger equation 

ih^\e{t),s) = H elec {t)\e{t),s), (11) 

with initial conditions \e(t = 0), s) = \e,s), the reduced density matrix automatically satisfies 
the correct equation of motion [Eq. ((HI)]. 

Within this framework, the equations for the nuclear trajectories are: 

• /.N Pn(t) 

M (12) 
p n (t) — — K (2u n (t) - u n+1 (t) - u n -i(t)) + 2aRe {p n>n+ i(t) - p n , n -i{t)} ■ 

The chain is taken to be clamped so that u\(t) = u^it) = and pi(t) = PN{t) = for all 
time, and Eq. (Tl2|) is valid for n = 2, • • • , N — 1. In turn, the orbitals that form p nm (t) 
satisfy Eq. (fTTl) . so that 



ih^-(n\e(t)) = [-t + a(u n+1 (t) - u n (t))} (n + l\e(t)) 

a/ " (13) 

+ [-t + a(u n (t) - u n ^(t))} (n - l|e(t)) 



for n, e = 1, • - • , N. Since the electrons are confined within the chain, (n|e(t)) = for 
n ^ {1, • • • , iV}. Equations (!T2l) and ( fTBl constitute a closed set of N(N + 2) coupled first- 
order differential equations that are integrated using an eighth-order Runge-Kutta method. 
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C. Nuclear initial conditions 



For the purpose of determining the nuclear initial conditions, the electronic state \<p{0)) = 
\E ) (H e [ ec \E ) = E \E )) is assumed to be well described by a single Slater determinant for 
which [recall Eq. 

M°)l4/v>(°)> = WM> ( 14 ) 

where f(e,s) is the initial electronic distribution (ground or excited) that takes values or 
1 depending on the initial occupation of each level with energy e and spin s. The starting 
optimal (minimum energy) geometry is obtained by minimizing the total energy of the chain 
by an iterative self-consistent procedure. Specifically, the energy gradient of the oligomer is 
given by 

dE(u) dH 

~du = ^( )Ift7~I^( )) = 2 « Re {Pm,m-l - Pm,m+l} + K(2u m - U m _i - U m+1 ) . (15) 

At the optimal geometry, for which the gradient equals zero, the m — 2, . . . , N — 1 displace- 
ment satisfies 

1 a 

U m = - (Wm+l + «m-l) - T7Re{p mjjn _i - p m ,m+l} ■ (16) 
L K 

Equation (fl6|) is solved iteratively with the additional constraint that the boundaries of the 
chain are clamped {u\ — Ujv = 0). 

Subsequently, a harmonic approximation to the nuclear ground-state wavefunction is 
obtained by performing a normal mode analysis around the equilibrium minimum energy 
geometry u° = (v,®, ■ ■ ■ ,u° N ) in the (ground or excited) initial electronic state \Eq). For 
this, the Hamiltonian is expressed as a sum of the static equilibrium configuration H and 
a dynamical part due to deviations from equilibrium 

H = H + H ' h + H' h , (17) 



where 



7V-1 



n=l,s 

N 2 K N ~ 1 
H U = E + Y E [ 2 ^ + l - U n)^ + 1 - Tin) + ( Vn+1 - Vn f 



n=l n=l 

with 7] n being the displacement of the n th monomer from its equilibrium position i] n = u n —u^. 
In order to get the potential energy of the chain around the equilibrium geometry, the 



quantity b is considered as a perturbation to H [40|, |4l| and we have to second order 
that, 



En — Ej 

(18) 

K y * 

+ y ^2l 2 ( U n+l - U° n )(Vn+l ~ Vn) + {Vn+1 ~ Vnf] 
n=l 

where we have traced over the electronic coordinates and assumed that the system is initially 
prepared in the electronic state \Eq) with energy Eq. Here {|v?i), E{\ are the eigenstates and 
eigenvalues of the A/"-particle electronic Hamiltonian in the optimal geometry Hq. 

A harmonic version of Eq. (fl8|) is obtained by making a Taylor expansion of the poten- 
tial around the equilibrium position and keeping terms up to second order in the nuclear 
displacements. We note that second-order perturbation in H^-ph is consistent with the har- 
monic approximation. The effective harmonic phonon potential energy thus obtained is: 



£ ham (T7) = E + l - Vnfn,mVm- (19a) 

n,m=2 

Here f n>m is the Hessian of the potential energy given by: 

d 2 E 



fn;, 



dr]n dr] % 



— Vnm + K{25 n ^ m — 5„ jm +l — 5 n ,m-l) ; (19b) 

r)=0 



where 



(19c) 



V 



n (e, e') = (e\n) ((n - l|e'} - (n + l|e')) + (n\e') ((e|n - 1) - (e|n + 1) 



In deriving Eq. (fl9l) we have imposed clamped ends on the polymer chain (rji — t]n — 0). 
The orbitals |e) and their associated energies e are obtained by diagonalizing the electronic 
Hamiltonian at the equilibrium geometry. The normal mode coordinates and frequencies are 



then computed by the standard analysis 42j. The eigenvectors of f nm provide the normal 



mode coordinates Qj(rj) and the associated eigenvalues Xj the normal mode frequencies 

U3j = y/Xj/M. 

A phase-space like description of the resulting nuclear quantum state is obtained by con- 
structing the associated nuclear Wigner phase-space distribution function p w (n, p). In the 



S 



normal- mode coordinates, pw(u,p) is just the product of the Wigner distributions associ- 
ated with each vibrational mode 

N-2 

Pw(u,p) = J] p 3 (Q 3 (u),P 3 (p)), (20) 

where Qj{u) is the normal mode coordinate of the j-th mode and Pj{p) its conjugate 
momentum. We take the chain to be initially prepared in its ground vibrational state so 
that |39j 

Pi{Qv p i) = \ exp(-Mco 3 Q 2 Jh) expi-Pf/frw.M) (21) 
Tin 

for j = 1, • • • , N—2. The 2N — 4 dimensional phase-space distribution in Eq. ( 1201) completely 
characterizes the initial quantum state of the nuclei. 

The ensemble of lattice initial conditions, {u l (0), p*(0)}, for the quantum-classical dynam- 
ics is obtained from a Monte Carlo sampling of the nuclear Wigner phase space distribution 
of Eq. (1201) . The average classical energy of the resulting ensemble coincides numerically 
with the zero-point energy of the lattice. The associated initial values for the orbitals 
{|e 1 )} are obtained by diagonalizing H e \ ec in the initial lattice geometries {u 1 }. Each initial 
condition i, together with the equations of motion, defines a quantum-classical trajectory 
(■u*(0),p*(0), \(p l (0))) — > (u l (t),p l (t), \(p l (t))) and the set is employed to obtain ensemble 
averages. Results shown here are averages over 10000 trajectories. 



III. RESULTS AND DISCUSSION 



Throughout, we study neutral oligomers with an even number N of CH units. In the 
ground state, the geometry of the chain consists of a centrosymmetric structure with perfect 
alternation of double and single bonds. The single-particle spectrum of chains of different 
length is shown in Fig. [2j It has a total width of 4t = 10 eV and consists of N/2 fully 
occupied valence states and N/2 initially empty conduction states. Note how the single- 
particle spectrum gets more dense as the number of CH units is increased. 

As a measure of electronic coherence and decoherence we follow the dynamics of the 
molecular polarization, defined by (/*(£)) = (\I/(t)|/t|\I / (t)), where denotes the vibronic 
wavefunction and where the dipole operator fi = ft e + fi^ has both an electronic fi e and a 
nuclear /In component. Doing so provides a measure of electronic coherence that is directly 
related to measurables. 
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FIG. 2: Orbital energies for optimized PA chains of varying length. The valence (conduction) 
orbital energies are plotted in blue (red). 



It is advantageous to connect this discussion of decoherence based on the polarization to 
standard measurements of decoherence Q, The density matrix associated with a general 
entangled vibronic Born-Oppenheimer state of the form = e~ zEnt ' h \(p n ) \Xn{t)} is 

given by 

Mt))m)\ = J2 e ^ nmt \^)\xn(t))^ m \( Xm (t)\, (22) 

nm 

where \ip n ) are the electronic eigenstates [H e i ec \ip n ) = E n \(p n )], \Xn(t)) the nuclear wavepacket 
associated with each electronic level and u nm = (E n — E m )/h. If our interest is in the 
electronic degrees of freedom only, then the vibrations are regarded as the environment. 
Since we have no interest in the behavior of the environment we trace over these modes to 
give the density matrix of the electronic subsystem: 



Pe (t) = J2 e ~ i0Jnmt (Xm(t)\Xn(t))\fn)(fr 



(23) 



Note that the off-diagonal elements of p e {t) are determined by the nuclear overlaps S nm (t) = 
(Xm(t)\xn(t)) and the loss of such coherences in p e {t) is a result of the evolution of the S nm (t) 
due to the vibronic dynamics. Standard measures of decoherence capture precisely this. For 
example, the purity of such entangled vibronic state is given by 



Tr(^(t)) = ^|(x m (t)IXn(t))P 



(24) 



and decays with the overlaps of the nuclear wavepackets in the different electronic surfaces. 

The polarization is also a useful measure of decoherence because its magnitude also 
depends on the S nm (t). To see this consider the expression for the polarization for the 



10 



entangled state in Eq. (122|) : 

(/*(*)) = E^(*)i^l^( t )) + E e "" nm Vr(Xm(t)iXn(t)}, (25) 

where /z™ n = (^ml/M^n)- Suppose that the PA chain is prepared in a spatially symmet- 
ric state where the initial nuclear state is invariant under reflection, i.e. pw(—u, —p) = 
pw(u,p). Since there is no symmetry breaking term in the Hamiltonian, this initial symme- 



try is maintained throughout the dynamics |44( and (Xn\^N\Xn) = AC™ = f° r an n - Under 
such conditions, the polarization 

(£(*)) = £ e--«™ Vr(Xm(*)|Xn(0> (26) 

is a direct measure of the off-diagonal matrix elements of the electronic reduced density ma- 
trix [cf. Eq. ([23])]. Its evolution and decay directly offers information about the decoherence 
dynamics. 

Note that in writing Eq. (I26p we have adopted the Franck-Condon approximation where 
the electronic transition dipole surfaces fi™ n (u) are assumed to depend weakly on the nuclear 
displacements. However even when this approximation is not valid, a decay in (/}(£)) will still 
signal a decay in the nuclear overlaps, albeit modulated by the dependence of the electronic 
transition dipoles on the nuclear coordinates. 

Thus, both Tr(pg(t)) and (fl{t)) are useful measures of decoherence and both decay with 
the overlaps of the nuclear wavepackets in different electronic states. The advantage of the 
polarization over the purity is that it is a physically accessible observable. Its limitation, 
however, is that (/*(£)) only signals coherences for which n™ n ^ 0. So, for instance, coher- 
ences between eigenstates of the same parity are absent in the polarization even when they 
would contribute to the purity. 

In the quantum-classical picture of the dynamics, the polarization is computed as an av- 
erage of the polarizations recorded for each of the Ai individual trajectories in the ensemble: 

I | M N 
i=l n=l 

where x l n (t) = (na + u l n (t)) is the position of site n at time t in the ith trajectory and e 
is the electron charge. The first term in Eq. ( 127|) comes from the dipole due to the nuclei, 
while the second one quantifies the electronic contributions. 
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There are two possible effects that can lead to a decay of the electronic coherences (re- 
call Fig. [1]): anharmonicities in the potential and population transfer into other electronic 
states. More precisely, if there is no population transfer into other states, and the electronic 
potential energy surfaces are bounded, then anharmonicities in the electronic potential en- 
ergy surfaces can lead to a spread of the nuclear wavepackets during evolution and thus to 
a decay of the nuclear overlap integral = (Xj{t)\Xi{t)) (wavepacket evolution in purely 
harmonic potentials lead to periodic recurrences in Sij and thus cannot lead to decoher- 
ence). Alternatively, population transfer into other electronic states can lead to decoherence 
by transferring population to states for which only poor overlaps of the nuclear wavepack- 
ets are possible with the states already involved in the superposition. This poor overlap 
arises because different electronic potential energy surfaces typically have substantially dif- 
ferent gradients and position of their minima in conformational space, leading to diverging 
evolution of the nuclear wavepackets in the excited state manifold. 

It should be noted that, traditionally, studies of the decoherence of a superposition state 
would not typically include loss of population from the state, which would be regarded as a 
relaxation, rather than decoherence, process. However, this distinction is meaningful when 
energy transfer and decoherence time scales are substantially different, the latter occurring 
on much shorter time scales than the former. Here, however, as shown below, population 
transfer between states occur rather quickly, making this subdivision less meaningful, and 
making such contributions quite significant in the time evolution of (/*(£)) and the purity. 
Here then, we use the term "decoherence" to relate to any process that causes loss of the 
coherence of (fi(t)). 

Below we discuss several examples of decoherence dynamics in PA chains. The interpre- 
tation of the results will be done in a wavepacket language and with the wavepacket picture 
of Fig. [I] in mind, even though the computations are performed in a mixed quantum-classical 
setting. Such wavepacket evolution is captured by the quantum-classical dynamics through 
the time dependence of the orbital energies and populations in the ensemble of trajectories. 
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FIG. 3: Electronic decoherence dynamics in chains of different lengths N. The figure shows the evo- 
lution and decay of the chain polarization when the system is initially prepared in a superposition 
between the ground and first excited state of the form in Eq. ([28]) . 



A. Decoherence between the ground and first excited state for chains of different 
lengths 

Consider first the decoherence dynamics of PA chains initially in a separable superposition 
state of the form: 

|*(0)> = -^(|^o) + ® |Xoo>, (28) 

where \(fo) is the ground electronic state, \<pi) is the first excited state (obtained by pro- 
moting an electron from the HOMO to the LUMO) and |xoo) is the ground state nuclear 
wavefunction in the ground electronic surface. Physically, such a superposition can be cre- 
ated by instantaneous (delta pulse) excitation of the relaxed ground state chain. Figure [3] 
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FIG. 4: Population of the LUMO orbital during the decoherence dynamics of a N-site SSH chain 
with the initial conditions of Eq. (|28p . 

shows (/!(£)) for chains with varying number of CH units (N). The high frequency oscil- 
lations in (fi(t)) are due to the difference in energy between the two states involved in the 
superposition (in this case the energy gap). The remaining time dependence arises from 
the wavepacket evolution in the excited state potential energy surface. For the four site 
chain, the polarization displays a fast initial decay with recurrences every ~ 30 fs. These 
recurrences arise from the time dependence of the overlap of the nuclear wavefunctions in 
the ground and excited electronic states [see Eq. (126]) ]. and signal the oscillatory motion 
of the nuclear wavepacket in the excited state potential. Between consecutive recurrences 
the amplitude of the polarization diminishes and eventually dies out, yielding a decoherence 
timescale of ~ 250 fs for N = 4. For N = 20 we observe only two of these recurrences, oc- 
curring every ~ 46 fs to yield a decoherence time of ~ 100 fs. For longer oligomers (N = 50 
and iV = 100) no recurrences are observed and the decoherence occurs in less than 10 fs. 

Additional insights into the decoherence dynamics can be extracted by considering the 
evolution of the population of the LUMO of the chain (Fig. H]). If the main decoherence 
mechanism is the anharmonic evolution of the nuclear wavepacket in the first excited state 
potential energy surface, then one should expect little population exchange with other levels 
of the chain. As shown in Fig. H] for N = 4 and N = 20 an almost negligible amount of 
population is transferred to other electronic states, suggesting that anharmonicities are the 
main source of decoherence. By contrast, for long chains (N = 50 and N = 100) the elec- 
tronic spectrum is so dense that a substantial amount of population is transferred from the 
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initially populated LUMO to other electronic states. This suggests that both anharmonic- 
ities and population decay to other electronic states contribute to the decoherence, leading 
to an evolution with no apparent recurrences. 

B. Decoherence of superpositions between excited states 

We now investigate how the decoherence dynamics changes when the initial superposition 
is between two excited states rather than between an excited and ground electronic state. 
For this we consider the two classes of model initial superpositions schematically represented 
in Fig. [5j In the first class, the initial state is of the form 

|*(O)) = i=(|^) + |0 i+1 ))®|xoo), (29) 

where = c\ s cjv/2,s|^o) (* £ {N/2 + 1, ■ • • , N}) is an electronically excited state obtained 
by promoting an electron from the HOMO to the zth orbital level of the ground state \tpo). In 
this superposition the initial nuclear state is taken to be the ground vibrational state in the 
ground electronic surface, |xoo)- Physically, such a superposition will arise via instantaneous 
excitation of the ground vibronic state into states \4>i) and \4>i+i), as depicted in the left panel 
of Fig. [51 By contrast, in the second class of superpositions the nuclei are taken to be initially 
prepared in the ground state distribution of the excited electronic state \<f>i), so that 

|*(O)> = -^(|0i) + |^i))®|xoi)- (30) 

The wavefunction \xoi) is obtained by finding the optimal geometry of the electronically 
excited state and then performing a normal mode analysis around this geometry, as discussed 
in Sec. Ill Cl Physically such a superposition will arise from instantaneous excitation of a chain 
vibrationally relaxed in state \<pi) to state \4>i+i), as represented in the right panel of Fig. |5l 

1. The case of a 20-site chain 

Consider first the case of a chain with 20 CH units. The left panels in Fig. [6] show the time 
dependence of the polarization when the system is initially prepared in the superposition 
in Eq. (|29l) for different z's. The figure shows that the decoherence dynamics can change 
substantially depending on the pair of states that are selected to form the superposition. 
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It is even possible to find superpositions for which the coherences are unusually long lived. 
For example, for i — 15 the electronic coherences survive for ~ 200 fs, a timescale that is 
comparable to the coherence lifetime observed in photosynthetic systems. 

Additional insights into the decoherence mechanisms are provided by the shape of (fi(t)) 
and by the dynamics of population in the excited orbitals (Fig. [7j). The polarization indicates 
that for % = 11 there is vibronic evolution in the excited states that leads to a decay and to 
the recurrence of the nuclear overlap integrals determining (/t(t)). By contrast, for i — 19 
this motion is not apparent in (//(£)) which shows a decay in ~ 60 fs with no apparent 
additional structure. The population dynamics (Fig. [7J upper panel) complements this 
picture by showing that for i — 11 a negligible amount of population decays to other levels, 
while for i = 19 the transfer of population to other levels is substantial. These observations 
suggest that for i — 11 the main mechanism for decoherence is due to anharmonicities in the 
excited state potential energy surfaces, while for i = 19 the main decoherence mechanism is 
due to population transfer to other electronically excited states for which only poor nuclear 
overlaps are possible. The case of i = 15 is discussed below. 

As an additional test of these observations consider the dynamics of superpositions be- 
tween the same set of levels but starting from Eq. (1301) . that uses a different initial nuclear 
state. The results are shown in the right panel of Fig. [6j As can be seen, for i — 11 chang- 
ing the initial nuclear state triples the coherence lifetime of the superposition, with three 
visible recurrences instead of one. Since for this superposition there is negligible amount of 
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*(0)) = -k(|&) + |&fi>)®lxo* 




100 150 200 250 300 



FIG. 6: Electronic decoherence dynamics in a 20 site chain for different initial superposition between 
excited states. The panels show the dynamics and decay of the chain polarization starting from a 
superposition of the form in Eq. ([29|) (left panels) or Eq. (|30|) (right panels). 

population being transferred to other electronic states (see Fig. [7]), the data confirms that 
anharmonicities in the potential energy surfaces of the excited states are the main source 
of the decoherence in this case. By contrast, for i — 19 changing the nature of the initial 
nuclear state has little effect on the decoherence dynamics, suggesting that the main deco- 
herence mechanism in this case is due to population transfer to other electronic states, as 
seen in Fig. [7J 

The case for i — 15 where long coherences are observed is different. For this superposition, 
little population is transferred to other electronic states and a change in the initial nuclear 
state has little effect on the decoherence dynamics. This suggests that this superposition 
is protected from decoherence both by the fact that the density of states is such that the 
two states involved in the superpositions are weakly coupled to other electronic states, and 
because the sampled potential energy surfaces are less anharmonic than in the other cases 
considered. 
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FIG. 7: Combined population of the i th and (i + l)*' 1 orbitals during the decoherence dynamics of 
a 20-site SSH chain starting from the initial superpositions depicted in Fig. [5j 




2. The case of a 100-site chain 

For larger systems the situation is qualitatively different. Figure |8] shows the dynamics 
of the polarization for chains initially in a superposition state of the form in Eq. (|29|) and 
Eq. fl30|) for different i. Figure |9] shows the associated change in population of the ith 
and (i + l)th orbitals. The electronic spectrum is so dense that upon evolution significant 
population is transferred to other electronic states. Irrespective of the type of superpositions 
considered or the initial nuclear state coherence decay in this chain is extremely fast, of the 
order of 50 fs. The electronic spectrum in this system is simply too dense to maintain 
electronic coherence. In (fi(t)), however, for i — 51 population loss is less than in the other 
two cases, consistent with the fact that some oscillatory character is visible for i = 51 in 
Fig. E 
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l*(o)> = + \<h+i)) ® Ixoo) |*(o)> = M^i) + \M) ® Ixoi) 




-1 1 1 1 1 1 1 1 -1 1 • 1 1 ■ 1 1 

O 50 100 150 200 250 300 50 100 150 200 250 300 

t (is) t (fs) 

FIG. 8: Electronic decoherence dynamics in a 100 site chain for different initial superposition states. 
The panels show the dynamics and decay of the chain polarization starting from a superposition 
of the form in Eq. (|29j) (left panels) or Eq. (I30j) (right panels). 

IV. CONCLUSIONS 

In this paper, we have presented numerical simulations of the electronic coherence dy- 
namics of PA oligomers of varying length in which the evolution of both electrons and nuclei 
are followed explicitly in a mixed quantum-classical approximation. We investigated the 
decoherence of superpositions, as manifest in the dynamics of the polarization, between the 
ground and excited and between pairs of excited states. Decoherence is caused by the decay 
of the overlap of the nuclear wavepackets associated with all electronic states involved in the 
superposition. Two basic mechanisms for such decay were identified: population transfer 
into other electronic states where only poor overlaps are possible, and vibronic evolution in 
anharmonic potential energy surfaces that lead to wavepacket spread. 

The simulations indicate that for long chains (e.g. N = 100) the electronic spectrum 
is so dense that decoherence is dominated by population decay into other states. In this 
case, no recurrences are observed in the polarization and the decoherence occurs in tens of 
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FIG. 9: Combined population of the ith and (i + l)-th orbitals during the decoherence dynamics 
of a 100-site SSH chain starting from the initial superpositions described in Fig. [5j 

femtoseconds. Further, the decoherence dynamics was found to be largely independent of 
the type of initial superposition that is subject to the decoherence. By contrast, for shorter 
chains (e.g. N = 20) the simulations indicate that the decoherence dynamics depends 
strongly on the initial vibronic state. We identified superpositions for which anharmonicities 
were the main source of decoherence and superpositions for which population transfer to 
other electronic states was determinant. Interestingly, we also observed a superposition 
state between excited states with coherence properties that are long lived, for ~200 fs. 
Such a superposition was found to be long lived because it is spectrally isolated from other 
electronic states and because the vibronic dynamics leads to a relatively slow spread of the 
nuclear wavepackets. 
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